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褐飞虱 发 生 的 相 空 间 线 性 回归 预测 模型 
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摘要 : 基于 混沌 理论 ， 通 过 对 江苏 省 太湖 地 区 农 科 所 褐飞虱 发 生 时 间 序 列 资料 的 分 析 ， 组 建 了 褐 飞 乔 Nilaparvata lugens Si 
发 生 时 间 序 列 的 相 空间 线性 回归 预测 模型 。 结 果 表 明 ，7 次 预报 的 平均 相对 误差 为 22.91%， 转 化 为 发 生 等 级 ， 预 测 准 确 率 
为 100 和 % ， 这 为 害虫 的 长 期 可 预测 性 提供 了 一 种 有 效 的 新 方法 。 
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Abstract: A regression model] of phase space for predicting brown planthopper (BPH) occurrence in the rice growing ar- 
ea of Taihu District, Jiangsu Province, was established on the basis of chaos theory. Its accuracy was very high with a 


relative error of only 22.91% . Therefore, the model provides a new and more efficient method for the prediction of the 


long-term occurrence of this insect pest. 
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福 飞 各 Nilaparvata lugens Stil 是 东南 亚 各 国 、 
日 本 、 朝 鲜 、 我 国 长 江 流域 及 以 南 稻 区 主要 迁 飞 性 
害虫 ， 近 二 十 多 年 来 ， 发 生 频 次 显著 增加 ， 危 害 程 
度 日 趋 严重 。 特 别 是 1987 年 褐飞虱 在 全 国 大 范围 
内 暴发 成 灾 ， 发 生 面积 占 全 国 水 稻 面 积 的 45%， 
损失 稻谷 200 万 吨 ，1991 年 ， 福 飞 虱 个 仅 在 全 国 南 
方 稻 区 大 发 生 ， 天 津 、 河 北 及 西南 部 分 省 份 的 稻田 
也 突 草 危害， 直接 经 济 损失 50 亿 元 人 民 币 ，1997 
年 又 遭受 全 国 性 的 宰 飞 各 大 发 生 ， 给 国民 经 济 和 农 
民 带 来 了 已 大 损失 〔 胡 伯 海 和 姜 瑞 中 ，1997)。 显 
然 ， 如 能 对 福 飞 者 的 灾害 发 生 作出 准确 及 超前 预 
报 ， 势 必 减 少 损失 ， 推 动 农业 可 持续 发 展 ， 造 福 人 
类 。 

统计 分 析 法 是 目前 褐飞虱 长 期 预测 工作 中 使 用 
较为 广泛 的 预报 方法 之 一 。 该 方法 计算 简单 、 快 
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速 ， 主 要 缺点 是 ， 当 预报 模型 不 合理 或 预报 因子 选 
取 不 当时 ， 预 报效 果 不 理 想 。 近 期 一 些 研究 论文 
(Ma et al.，2000; 马 飞 等 ，2001) 不 仅 报道 了 袜 飞 
恒 发 生 系 统 具 有 分 形 性 质 和 混沌 特性 ， 而 且 还 较为 
系统 地 阐述 了 褐 飞 重 发 生 系统 的 可 预报 尺度 分析 
了 长 期 预测 准确 性 不 稳定 的 内 在 原因 。 随 着 混沌 理 
论 和 应 用 技术 研究 的 不 断 深入 ,混沌 时 间 序 列 的 建 
模 和 预测 已 成 为 混沌 信息 处 理 领域 中 的 研究 热点 
《 王 明 进 和 程 乾 生 ，199%7; 程 瑶 生 等 ，2000; 张 家 
树 和 肖 先 赐 ，2000)。 作 者 基于 相 空 间 重 构 理 论 ， 
提出 了 用 于 时 间 序 列 预测 的 相 空 间 线性 回归 预测 模 
型 ， 并 用 江苏 省 太湖 地 区 农 科 所 调查 记载 的 神 飞 乱 
发 生 量 数据 资料 对 该 模型 的 预测 效果 进行 了 检验 。 
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许 晓 风 等 .褐飞虱 发 生 的 相 空间 线性 回归 预测 模型 


1 原理 与 方法 


1.1 相 空 间 重 构 

设 时 间 序 列 人 (5 )}， 其 时 间 间 隔 是 At (单位 
时 间 )。 显 然 ， 吸 引子 的 结构 特性 就 包含 在 这 时 间 
序列 中 。 为 了 从 实际 的 单 变 量 时 间 序 列 中 估计 出 系 


xCt1) xC ts) 
xlt1+t) wk 年 这 2 
xCt) + 27) 了 所 + 27) 


CE 


XH + Cd - Dr) 
Yl) 


和 有 + Cd — 1)r) 
Y(z,) 


(1) 式 的 各 列表 示 了 4 维 空 间 的 各 相 点 Y(z,)《 见 
(2) 式 ) 的 & 个 分 量 值 ， 分 别称 为 第 1，2，… 及 第 
d 个 分 量 。 式 中 的 z 为 延 滞 时 间 ， 它 是 tx(z,) 序 
列 的 间隔 Az 的 整数 倍 ，d 为 整数 。 通 常 d =3，4， 
5$，6，…。r，d 的 取 值 可 在 模式 用 于 不 同 参数 预 
报 试 算 中 确定 。 相 点 数 m=mrm- (d-1) r。 

做 预报 ， 即 是 想 知道 下 一 时 刻 的 xCz, +1)》 值 ， 
由 (1)，(2) 可 知 ，Y(Czt, + 1) 的 4 个 分 量 中 只 
x(Ct, +1)》 是 未 知 的 ， 它 正 是 我 们 待 求 的 量 。 
1.2 ” 相 空 间 线 性 回归 模式 的 建立 

某 一 参考 态 Y(s,〉 附近 有 许多 相 轨 道 ， 找 出 
Y (1 ) 的 几 个 最 近邻 态 Y,,(z1,) 并 考查 其 所 处 相 轨 
道 的 统计 行为 《关系 )， 利 用 相似 原理 将 其 移植 于 
参考 态 Y(1, ) 所 处 相 轨 道 ， 以 求 得 下 一 时 刻 的 参量 





b 
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统 的 信息 ，Packard〈1980》 提出 了 用 时 间 序 列 重 构 
吸引 子 的 相 空 间 图 和 象 。 
Takens〈1981) 提出 稻 入 定理 ， 表 明 只 要 相 衬 
间 的 维 数 d 足够 大 〈 通 常 为 2D,+ 1)， 它 就 可 以 刻 
划 D。 维 混沌 吸引 子 。 
把 时 间 序 列 {x《i,)} 延 拓 成 gd 维 相 空 间 的 一 
个 相 型 分 布 ; 


wk Bd dz) 
2 ‘xlt, -Cd 2)7) 
ek A ‘xlt, -Cd -30r) (1) 
xt + Cd Lr) X(E 7 
“Cp ) CC (2) 
x《t, +1)。 具 体 过 程 如 下 : 
首先 找 出 参考 Y(t: ) 的 最 近邻 态 了 (5 ) 
Yatt) = min| Yl#,) — Ys) (3) 


j= 1,2,A(m - 1) 
式 中 min 表示 最 小 ，| | 表 欧 氏 模 。 由 (3》 可 以 
找到 最 小 、 次 小 …… 个 近邻 点 Y(t, )。 然 后 考查 
这 个 相 轨 道 的 统计 行为 。 了 Y,(#,) 经 时 间 TCT = 
At，2AtsA》 后 ， 在 其 轨道 上 演变 成 Y, (zs, + 了 )。 
在 7 不 太 大 的 情况 下 《 即 相 轨道 演变 极 短 )， 假 设 
Yi 《t+ Ty 和 Y(t,) 服 从 以 下 的 线性 关系 ; 
| Yi + Th) = a + 6 YCs), | 
用 最 小 二 乘法 可 以 得 到 
> ,| Ys, + TD 
SS 天 


人 


(4) 





- b>, | YC) (5) 





kD | 全 | | kd + 下 | - 2 | Ye | A + TP (6) 
. ky Ye CDs | Ys | 
从 记 个 YCz) 及 上 个 Y, 《zt, + T) 的 有 关 资 料 里 (7) 式 右 端 之 a、56 由 (5)、(6) 给 出 ，Y 了 Y(t, ) 由 


[它们 完全 由 (1)、(2) 确定 ]， 求 得 系数 a、&b。 
将 近邻 相 扫 道 的 相似 行为 《4) 移植 于 参考 态 
Y 《zt,》 所 处 的 相 轨 道 ， 则 预报 方程 为 
| Yt, + T) | a+b| Ys,)| (7) 
在 相 空间 里 ， 我 们 可 以 逐一 来 考查 不 同 相 点 在 
不 同时 刻 的 演化 情景 ， 使 预报 模式 的 函数 形式 或 系 
数 可 因 时 而 变动 。(7) 式 给 出 的 线性 形式 在 各 相 点 
虽 不 变 , 但 系数 a、6 却 随 相 点 移动 而 变化 ， 这 就 
避免 了 传统 统计 模式 的 缺陷 。 


(1)、(2) 给 定 ， (7) 式 左 端的 YC + T) 中 具有 x 
(5 + T) 分 量 是 未 知 的 ， 由 《7) 式 可 求 出 x( + 
7)。 了 TT 是 序列 tx (zi,) 的 时 间 间 隔 Ar 的 整数 倍 ， 
它 可 随 预 报 要 求 而 定 。 


2 方法 的 应 用 


研究 中 使 用 的 神 飞 筷 发 生 资料 ， 取 自 种 植 单 季 
晚稻 为 主 的 江苏 省 太湖 地 区 农 科 所 〈120.30 EE， 
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31.50° N》1986 ~ 1998 年 13 年 间 7~9 月 田间 观察 
男 系 统 调查 资料 ， 序 列 长 度 为 381。 在 进行 混沌 时 
间 序 列 相 空间 线性 回归 预测 模型 的 建立 时 ， 我 们 采 
取 了 以 年 为 间隔 抽取 、 以 旬 为 单位 等 间隔 数据 点 的 
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方法 。 以 积累 的 7 月 下 旬 的 数据 可 预测 下 一 年 7 月 
下 旬 的 袜 飞 虱 发 生 量 ， 同 样 以 积累 的 8 月 上 旬 的 数 
据 可 预测 下 一 年 8 月 上 旬 的 福 飞 剧 发 生 量 ， 依 次 类 
推 可 以 得 到 全 年 7~ 9 月 各 旬 的 褐飞虱 发 生 量 预测 。 


表 1 1998 年 7~9 月 各 和 旬 褐 飞 虱 发 生 量 的 预测 结果 检验 
Table 1 Test of the predictive accuracy of the model on N. Iugens (BPH) occurrence between July and September, 1998 


月 .x 日 实际 值 预测 值 相对 误差 绝对 值 (%》 
Month. Day Observed value Forecasted value Absolute value of relative error Simulated occurrence degree 
7.30 36 39.6735 10.20 
8.10 140 109.1692 22.02 
8.20 266 162.5923 38.88 
8.30 250 173.460 30.62 
9.10 1 630 1 860 14.11 
9.20 4 460 5 512 23.59 
9.30 5 340 6 458 20.94 
平均 误差 average error 22.91 


实际 发 生 程度 检验 
Observed occurrence degree Test 
轻 发 生 light occurrence ~ 

轻 发 生 light occurrence 
轻 发 生 light occurrence 
轻 发 生 light occurrence 


中 发 生 medium occurrence 


模拟 发 生 程度 * 


轻 发 生 light occurrence 
轻 发 生 light occurrence 
轻 发 生 light occurrence 
轻 发 生 light occurrence 
中 发 生 medium occurrence 


大 发 生 SEVEIE OCCLITENCE 大 发 生 SeEVEIe OCCUITENCE 


< 和 < < < 人 < 人 


大 发 生 SEVEIE 0OCCUHITEICE 大 发 生 SeEVEIe 0OCCUITECe 


* 根据 当地 褐飞虱 发 生 程度 等 级 划分 标准 百 穴 虫 量 < 500 头 ， 轻 发 生 ， 百 穴 虫 量 在 500 ~ 1 000 头 ， 偏 轻 发 生 ;， 百 穴 虫 量 在 1 000 ~ 2 000 
头 ， 中 等 发 生 ;， 百 穴 虫 量 在 2 000 ~ 3 000 头 ， 偏 重 发 生 ;， 百 穴 虫 量 >3 000 头 ， 大 发 生 

x* OQccurence degree of BPH was based on local classification standards: 1) BPH number/100 rice plants (B/100R) < 500, light occurrence; 2) B/100R = 500 
-1 000, a little more than light occurrence; 3) B/100R = 1 000 -2 000, medium occurrence; 4) B/100R = 2 000— 3 000, a little more than medium occur- 


rence; 5) B/I00R > 3 000，severe occurrence 


我 们 构造 了 10 维 相 空间 ， 通 过 计算 各 点 与 
Y(t ) 之 间 的 欧 氏 距离 找 出 了 (5 ) 的 最 近邻 域 点 大 
为 7 个 。 根据 公 式 (4)、(5)、(6)、(7) 得 到 1998 
年 7~9 月 各 旬 福 飞 重 的 发 生 量 ， 结 果 见 表 1。 从 
表 1 给 出 的 1998 年 7~9 月 各 旬 褐 飞 恒 发 生 量 的 实 
际 情况 与 预报 结果 的 对 比 中 不 难看 出 ， 预 报 结果 是 
令 人 满意 的 ，7 次 的 预报 平均 相对 误差 为 22.91% ， 
转化 成 的 发 生 程 度 与 实际 发 生 程度 相符 ， 预 报 准 确 
率 达 100%6 。 


3 结论 与 讨论 


(1)》 本 文 作者 基于 混沌 动力 系统 的 相 空间 延 返 
坐标 重 构 ， 利 用 混沌 序列 固有 的 确定 性 和 非 线 性 ， 
提出 了 用 于 混沌 时 间 序 列 预 测 的 相 空间 线性 回归 预 
测 模型 ， 通 过 对 1998 年 7~ 9 月 各 旬 的 褐飞虱 发 生 
量 实测 数据 资料 进行 预报 ， 得 到 了 7 次 预报 的 平均 
相对 误差 仅 为 22.91 和 多 ， 转 化 为 发 生 程度 与 实际 发 
生 程度 相符 ， 预 报 准确 率 达 100% 的 预报 结果 。 这 
就 为 害虫 时 间 序 列 的 可 预测 性 及 新 型 预测 模型 的 构 
建 提供 了 可 靠 的 依据 。 

(2》 相 空间 线性 回归 预测 方法 要 求 数据 愈 多 愈 


好 ， 本 试验 仅 用 了 1986 ~ 1998 年 的 历史 资料 ， 数 
据 总 数 不 足 400 个 ， 相 点 数 m = n -(d -1l)r 也 在 
400 以 内 。 故 相 轨 道 不 是 很 稠密 ， 影 响 近 邻 点 的 选 
择 以 及 a、&b 的 计算 ， 最 终 影响 预报 的 精度 。 如 果 
能 够 获得 更 多 的 资料 ， 用 于 实时 预报 ， 精 度 将 会 进 
一 步 提 高 ， 误 差 大 于 22.91% 的 份额 也 将 会 降低 。 

《3) 混沌 时 间 序 列 预 测 方法 是 一 新 的 很 有 发 展 
前 途 的 预测 方法 。 本 文 作者 提出 的 相 空 间 线性 回归 
预测 模型 ， 具 有 其 独特 的 优点 ， 主 要 表现 在 能 够 有 
效 地 减 小 测报 误差 。 在 实际 应 用 过 程 中 ， 随 着 对 害 
虫 发 生 系统 内 部 运行 机 制 的 进一步 揭示 ， 长 期 预报 
一 一 这 一 长 期 困扰 害虫 测报 的 难题 将 会 得 到 较 好 的 
解决 。 
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